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Abstract 



Perturbative corrections of the order of to inclusive double differen- 
tial lepton distribution from b quark decay are considered. A perturbative 
correction to the charged lepton energy spectrum has been calculated for an 
arbitrary charged lepton mass. The perturbative contribution suppresses the 
partial rate but almost does not change the shape of energy distribution. Ap- 
plications of our result to semileptonic B meson decays are briefly discussed. 
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I. INTRODUCTION 



The precise determination of weak mixing angles V^t and Vub is a demanding task. In 
spite of progress in this field they still remain ones of the worse known parameters of the 
Standard Model. The uncertainties which appear here are both of an experimental and the- 
oretical origin. The relatively large theoretical errors mainly reflect the lack of quantitative 
knowledge about the structure of hadrons and QCD higher order perturbative corrections 
to the amplitudes of weak decays of b quarks. The most valuable source of information 
about the weak mixing angles are the semileptonic decays of B and B mesons. The lep- 
tons in the final state do not interact strongly and the process is less affected by unknown 
QCD effects than a hadronic decay. Furthermore pseudoscalar B mesons are the simplest 
bottom hadrons. The b quark mass is about 5 GeV and thus it exceeds roughly ten times 
typical energy scales which characterize the infrared dynamics in the hadrons. Moreover the 
presence of this mass justifies the perturbative treatment of most of the processes involving 
the b quark. The simple facts have given rise to a quantitative description of dynamics of 
hadrons containing heavy quarks (Heavy Quark Effective Theory In the framework 

of HQET many observables describing heavy hadrons may be expressed as a power series in 
l/rrih. In particular it was shown in Ref. 0, that the inclusive lepton distributions from a 
bottom hadron decay may be treated in such a way. It follows from the operator product 
expansion (OPE) that a matrix elements which should be evaluated to derive the distribu- 
tions may be expanded into a series of local operators characterizing the decaying bound 
state. The very advantage of this approach is that subsequent unknown non-perturbative 
matrix elements are suppressed by increasing powers of m?,. The leading term corresponds 
to a parton contribution to the process. As argued in Ref. the next-to-leading term van- 
ishes. The I /ml corrections have been calculated by for a case of massless and [p!0|-[r2 
massive lepton. Recently also the third order terms have become known W^ . 



The first order perturbative QCD corrections to the inclusive lepton distributions in a 
process of decay: b qW are as important as the HQET corrections for the corresponding 
B decay. They have been evaluated [14,15| for the vanishing lepton mass. In the case of 
a non-zero lepton mass only a differential distribution of the lepton pair invariant mass is 



known to the first order in strong coupling constant |T^. In the present article we present 
our recent calculation of the first order QCD correction to the double differential inclusive 
lepton distribution from b decay with a massive lepton in the final state. The complete 



analytical result and details of the calculation will be published elsewhere |T3]. Here we 
give results for the perturbative correction to the r lepton energy spectrum which has been 
obtained by numerical integration of this double differential distribution. 



II. KINEMATICAL VARIABLES 

The purpose of this section is to define the kinematical variables which are used in this 
paper. We describe also the constraints imposed on these variables for three and four-body 
decays of the heavy quark. 

The calculation is performed in the rest frame of the decaying b quark. Since the first 
order perturbative QCD corrections to the inclusive process are taken into account, the final 
state can consist either of produced quark c, lepton r and r anti-neutrino or of the three 
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particles and a real gluon. The four-momenta of the particles are denoted in the following 
way: Q for b quark, q for c quark, r for the charged lepton, u for the corresponding anti- 
neutrino and G for the real gluon. By the assumption that all the particles are on-shell, the 
squares of their four-momenta are equal to the squares of masses: 

2 2 2 2 2 2 r~<2 n /i\ 

The four- vectors P — q + G and W — t + u characterize the quark-gluon system and the 
virtual W boson respectively. We define a set of variables scaled in the units of mass of 
heavy quark m^: 

_m| _'2Er _ _ 

ml ml ml m^ m^ 

We introduce hght-cone variables describing the charged lepton: 

T± = ^(x±^x2-47?) (3) 

The system of c quark and real gluon is characterized by the following quantities: 

Ps{z) = ,/^~~z=^[l + f + z'-2{t + z + tz)Y/\ 
P±{z) = Po{z)±Ps{z), 

(4) 

where Po{z) and Psiz) are the energy and length of the momentum vector of the system in 
b quark rest frame, yp{z) is the corresponding rapidity. Similarly for virtual W: 

Wo(z) = l{l + t-z), 
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W:iiz) = -t = -[l + t' + z^-2{t + z + tz)Y/\ 

W±{z) = Wo{z)±W^{z), 

(5) 

^From kinematical point of view the three body decay is a special case of the four body one 
with vanishing gluon four-momentum, what is equivalent to ^ = ^. It is convenient to use 
in this case the following variables: 
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Y,^y,(g)^l\n^, Y^^yM = \^^—- (6) 

We express also the scalar products which appear in the calculation by the variables x, t 
and z: 

Q-P^^{l + z-t) r-u=l{t-r]) 

Q-iy =1(1- z-x + t) T-P =\{x-t-r]) (7) 

Q-T^\x v-T = \{l - X - z + rj) 

All of the written above products are scaled in the units of the mass of h quark. 

The allowed ranges of x and t for the three-body decay are given by following inequalities: 

2^/?? <a:;<l + 77-^ = Xmax, (8) 

^. = r_{l-^)<t<r,{,-^yt. (9) 

(a region A). In the case of the four-body process the available region of the phase space is 
larger than the region A. The additional, specific for the four body decay area of the phase 
space is denoted as a region B. Its boundaries are given by the formulae: 

207 <X< Xmax, V <t <h (10) 

We remark, that if the charged lepton mass tends to zero than the region B vanishes. 

One can also parameterize the kinematical boundaries of x as functions oit. In this case 
we obtain for the region A: 

r7<^<(l-^)^ w_ + ^<x<w+ + ^, (11) 
and for the region B: 

?7<i< V^(l-Y^)' 2^<x<w^ + ^. (12) 

The upper limit of the mass squared of the c-quark — gluon system is in both regions 
given by 

Zmax = (1 - - t/T+), (13) 

whereas the lower limit depends on a region: 

{Q in the region A , . 

(1 - r_)(l - t/r_) in the region B. ^^^^ 
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III. EVALUATION OF THE QCD CORRECTIONS 

The QCD corrected differential rate for 6 — > c + r~ + reads: 

dr = dFo + ciri,3 + ciri,4, (15) 

where 

dTo = Glml\VcKM\^Mo,^dn3{Q;q,T,p)/7r' (16) 
in Born approximation, 

dr,,s = '^asGlml\VcKM\''Ml:,dns{Q;q,T,u)/7r'' (17) 

comes from the virtual gluon contribution and 

c?ri,4 = '^asGlml\VcKM\^M]:^^dn4iQ;q,T,iy)/Tr^ (18) 

describes a real gluon emission. Vckm is the Cabbibo-Kobayashi-Maskawa matrix element 
associated the 6 to c or m quark weak transition. Lorentz invariant n-body phase space is 
defined as 

d7^„(P;pl, . . . ,pn) = S^'\P - n (19) 

i * 

In Born approximation the rate for the decay into three body final state is proportional to 
the expression 

Mq^^^ Fo{x,t) =Aq-T Q-v ={1- g-x + t){x-t-ri), (20) 

where the quantities describing the W boson propagator are neglected. Interference between 
virtual gluon exchange and Born amplitude yields: 

■^1,3 = -[q-rQ-iy Ho + gQ-vQ-T H++ q-u q-r H_ + 
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+ -g v-r {H+ + H_) + -r)gQ-u {H+ ~ H_ + Hl) - -rj q-u H^], 



(21) 



where 

Ho = 4(l-ypPo/P3)lnAG + (2po/P3) 



LiJl-^) 



- Li2 (^1 - ^ j - + 1) + 2(ln ^ + Yp) {Y^ + Yp) 
+ [2psYp + {l-g-2t) In ^/t + 4, 
= ^[l±{l-g)/t%/ps±j]n^, 

Hl = ^(l-lnV^) + ^lnV^+|F,p3 + f^. 



(22) 



In M-i -^ infrared divergences are regularized by a small mass of gluon denoted by Ac- Ac- 
cording to Kinoshita-Lee-Naunberg theorem, the infrared divergent part should cancel with 
the infrared contribution of the four-body decay amplitude integrated over suitable part of 
the phase space. 

The rate from real gluon emission is proportional to 



M 



Bo 



1,4 



+ 



(g-G)2 Q-GP-G (P-G) 



(23) 



where 



B- = q-T [Q-u {Q-G - 1) + G-u - Q-uQ-G + G-u Q-G], 

B- = q-T[G-uQ-q - q-uQ-G + Q-u {q-G - Q-G -2q-Q)] + 

+ Q-V {Q-rq-G - G-rq-Q), 
B- = Q-u{G-Tq-G -QT-P). 



(24) 



Integrating and adding all the contributions one arrives at the following double differen- 
tial distribution of leptons: 



dV 



12ro 



dxdt \ l2To^F^{x,t) 



Fo{x,t) - ^Ff{x,t)\ for {x,t) in A, 
''^TrBf^+\ for(x,t)inB 



where 



" 1927r3 
Fq{x, t) = {1 — g — X + t){x — t — T]) , 



(25) 

(26) 
(27) 



and the functions F^{x,t), Ff(x,t) describe the perturbative correction in the regions A 
and B. Explicite formulae for F^(x,t) and F^(x,t) will be given in [|l^. The factor of 12 
in the formula ( pSj) is introduced to meet widely used convention for -^0(2^) and 

The obtained results were tested by comparison with earlier calculations. One of the 
cross checks was arranged by fixing the mass of the produced lepton to zero. Our results are 
in this limit algebraically identical with those for the massless charged lepton On 
the other hand one can numerically integrate the calculated double differential distribution 
over X, with the limits given by the kinematical boundaries: 



W++V/W+ dY dV . . 

-r^{x,t; g,r]) = —{t; g,r]) 
2^ dx dt dt 



(28) 



Obtained in such a way differential distribution of t agrees with recently published |T6 



analytical formula describing this distribution. This test is particularly stringent because 
one requires two functions of three variables {t, g and t]) to be numerically equal for any values 
of the arguments. We remark, that for higher values of t only the region A contributes to 
the integral ( PE| ) and for lower values of t both regions A and B contribute. This feature of 
the test is very helpful — the formulae for Fi{x, t), which are different for the regions A and 
B can be checked separately. 
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FIG. 1. (a) The distributions fo{x), fi{x) and (b) the ratio fi{x)/fo{x) for the pole mass of 
the b quark nib = 4.5 GeV (dotted), nib = 4.75 GeV (sohd) and nib = 5.0 GeV (dashed). 



IV. DIFFERENTIAL DISTRIBUTION OF r ENERGY 



The point of interest to check how the QCD corrections change energy spectrum of the 
charged lepton. This aim may be reached by integration of the double differential lepton 
distribution over the lepton pair invariant mass: 

dt, (29) 



dx J Tj diX dt 

where ^2 is the upper kinematical boundary for t given the formula (|TI|) . The decomposition 
of the resulting distribution into the Born term and the perturbative QCD correction yields 
in a natural way definitions of functions /o(a;) and /i(x): 



dr 

T- = i2ro 

dx 



(30) 



The analytical formula for fo{x) reads 



fo{x) = 2^x^ - Arj [xl[x^ - 3x(l + v) + H + 

+xl[-3x^ + 6x{l + ri) - 12ri]], 



(31) 

where following we introduced 

xo = l- g/il + v-x). (32) 
An equivalent expression for /o(x) is 



2 



fo{x) = ^xJx"^ - 4:1] I — ^1 [3(1 + 1]) -2x + g- Ari/x + 

\1 + 7] — X J 

+2q{1 + 7]- 4:i]/x)/{l + 1]- x)], 



(33) 



6 



where Xmax is given by (§). The latter formula clearly exhibits the behavior of /o(a;) for x 
close to the upper kinematical limit. 

The integration of F^'^{x,t) was performed numerically for different masses of 6-quark 
with fixed rrib — nic = 3.4 GeV and m^- = 1.777 GeV. The functions fo{x) and fi{x) for 
nib = 4.75 GeV are plotted on Fig. la and the ratios /i(x)//o(x) for three different realistic 
values of rrib are plotted on Fig. lb. As can be easily seen the ratios have logarithmic 
singularities at the upper end of the spectra. Such a behavior would lead to a inconsistence. 
The standard solution to problems of this kind is an exponentiation which yields well known 
Sudakov form factor [IE\. Far from the end point the ratio of the correction term to the 
leading one is almost constant and close to 2. It means that the perturbative correction 
changes rather the normalization than the shape of lepton energy distribution. 
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FIG. 2. The QCD corrected r lepton spectrum from the b quark decay for different values of 
as- The mass of b quark is chosen as 4.75 GeV. 



The obtained distributions of the scaled charged lepton energy for rUb = 4.75 GeV 
with and without perturbative QCD corrections are shown in Fig. 2. The strong coupling 
constant was chosen as 0.2 and 0.4 since the energy scale for this process in not known until 
the second order QCD corrections are evaluated. The value of for this decay is expected 
to lay between the two numbers. 

The knowledge of perturbative corrections to lepton energy is essential for fixing HQET 
parameters, especially Ai and A Especially analysis of moments of the lepton energy 

spectrum and other quantities involving integration over the energy distribution appeared 
particularly valuable for this purpose |jl9l as was earlier suggested in Refs. • 
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V. CONCLUSIONS 



The first order QCD corrections to the double differential inclusive lepton distributions 
from b quark semileptonic decay have been calculated for a massive fermion in the final 
state. Non-trivial cross checks of the final the result have been performed. We remark that 
including a real gluon radiation on the parton level yields a increase of the phase space 
available in the decay process. The QCD corrected r energy spectrum has been obtained. 
The effect of the correction may be estimated as about 10% of the magnitude of uncorrected 
distributions. 

The presented above results can be utilized to improve an analysis of semileptonic decays 
of beauty hadrons with a r in the final state. Thus the values of involved weak mixing 
angles may be fixed more exactly. The decrease of theoretical uncertainty increases the 
sensitivity to hypothetic deviations from the Standard Model [pl] -p3| which should have to 
be particularly distinct in the case of the heaviest family. The better understanding of the 
perturbative QCD effects allows one to perform more stringent tests of HQET predictions 
|T0HT2| and narrow the error bars for HQET parameters. Moreover one can extract more 
precisely some information about masses of quarks and strong coupling constant from the 
data. Finally, the process that we considered may appear a background for other processes 
so precise theoretical knowledge about the process is valuable. At present however the 
statistics of measured b c{u)ti'j. transitions is rather low and ten-percent effects are 
not seen. Probably the application of provided here formulae to the expected data from 
i?-factories will be really fruitful. 
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